arXivil 506.06649V 1 [nlin.CG] 22 Jun 2015 


An example of degenerate hyperbolicity in 
a cellular automaton with 3 states 


Henryk Fuks and Joel Midgley-Volpato 

Department of Mathematics and Statistics, Brock University 
St. Catharines, Ontario L2S 3A1, Canada 


Abstract. We show that a behaviour analogous to degenerate hyperbol¬ 
icity can occur in nearest-neighbour cellular automata (CA) with three 
states. We construct a 3-state rule by “lifting” elementary CA rule 140. 
Such “lifted” rule is equivalent to rule 140 when arguments are restricted 
to two symbols, otherwise it behaves as identity. We analyze the structure 
of multi-step preimages of 0, 1 and 2 under this rule by using minimal fi¬ 
nite state machines (FSM), and exploit regularities found in these FSM. 
This allows to construct explicit expressions for densities of Os and Is 
after n iterations of the rule starting from Bernoulli distribution. When 
the initial Bernoulli distribution is symmetric, the densities of all three 
symbols converge to their stationary values in linearly-exponential fash¬ 
ion, similarly as in finite-dimensional dynamical systems with hyperbolic 
fixed point with degenerate eigenvalues. 


1 Introduction 


In a linear continuous-time dynamical system given by x = Ax, if x : M —>■ M" 
and A is a real n x n matrix with all eigenvalues distinct and having negative real 
parts, x(t) tends to zero exponentially fast as t —)■ oo. The same phenomenon 
can be observed in nonlinear system x = f(x) (where f : M" —>■ M") in a vicinity 
of hyperbolic fixed point, as long as the Jacobian matrix of f evaluated at the 
fixed point has only distinct eigenvalues with negative real parts. If, on the other 
hand, the matrix A has degenerate (repeated) eigenvalues, the convergence to 
the fixed point can be polynomial-exponential, that is, of the form 
where P{t) is a polynomial and 5 > 0. 

In discrete-time dynamical systems things are quite similar. For example, the 
linear system system 
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is defined by a matrix which has degenerate (double) eigenvalue |, thus polynomial- 
exponential (linear-exponential in this case) convergence to the fixed point (0,0) 
is expected. Indeed, the solution is 
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and we can clearly see the aforementioned linear-exponential convergence. 

Cellular automata are infinitely-dimensional dynamical systems, yet a be¬ 
haviour similar to hyperbolicity in finite-dimensional systems has been observed 
in many of them. In particular, in some binary cellular automata in one dimen¬ 
sion, known as asymptotic emulators of identity^ if the initial configuration is 
drawn from a Bernoulli distribution, the expected proportion of ones (or zeros) 
tends to its stationary value exponentially fast [3]. 

Furthermore, an example of a probabilistic CA has been recently found [5] 
where the density of ones converges to its stationary value in a linear-exponential 
fashion, just like in the case of degenerate hyperbolic fixed points in finite¬ 
dimensional dynamical systems. Could such behaviour be observed in deter¬ 
ministic CA as well? The purpose of this paper is to provide an example of such 
deterministic CA. 

We will consider 3-state nearest-neighbour CA obtained from elementary 
binary CA by “lifting” them to 3-states. What we mean by this is the following 
construction. Let g : {0,1}^ —>■ {0,1} be a local function of elementary CA 
satisfying g{0,0,0) = 0 , g(l, 1 , 1 ) = 1 , and let fg : { 0 , 1 , 2 }^ —>• { 0 , 1 , 2 } be 
defined by 

{ g{xi,X2,X3), xi,a; 2 ,a :3 G (0,1} 

25 a;i,a; 2 ,a :3 G {0,2} 

g{xi -1,X2- 1, xa - 1) -I- 1, Xi,X 2 , 3:3 G {1, 2} 

X 2 , otherwise. 

This construction ensures that when fg is restricted to two symbols only, it 
becomes equivalent to g, otherwise it behaves as identity. Conditions 5 ( 0 ,0,0) = 
0 , (;(!, 1 , 1 ) = 1 ensure that there are no conflicts, so that, for example, /( 2 , 2 , 2 ) 
is the same no matter if we apply second or third case of eq. (|^ . 

We studied dynamics of fg for a number of elementary rules g. One of the 
most interesting of them is the case of g being elementary CA rule with Wolfram 
number 140, defined as 

g{Xi,X2jX3) = X2- XiX2-^XiX2X3, (4) 

where a;i,a; 2 ,a ;3 G {0,1}. As we will see, it actually exhibits degenerate hyper¬ 
bolicity. In what follows, we will refer to fg with g given by eq. Q as “rule 
140”. We will also drop the index g and refer to fg simply as /. An example of 
a spatio-temporal pattern produced by this rule is shown in Figure 

Let us first introduce the notion of density polynomials. Let A = {0,1,2}. A 
finite sequence of elements of A, b = 6162 ..., will be called a block (or word) 
of length n. The set of all blocks of elements of A of all possible lengths will be 
denoted by A*. 

A block evolution operator corresponding to / is a mapping f : A* 1 —>■ A* 
defined as follows. Let a = 0302 ... a„ G A" where n > 3. Then f(a) is a block 
of length n — 2 defined as 

f(a) = /(oi, 02 , a 3 )/(a 2 , 03 , 04 )... /(a„_ 2 , o„_i, a„). (5) 


fg{xi,X2,X3) = 


2 


i 



Fig. 1. Sample spatio-temporal pattern generated by 3-state rule 140. White, lighter 
gray and darker gray cells (blue in color version) correspond, respectively, to 0, 1 and 2. 


If f(b) = a, than we will say that b is a preimage of a, and write b G f“^(a). 
Similarly, if f"(b) = a, than we will say that b is an n-step preimage of a, and 
write b S f“"(a). 

Let the density polynomial associated with a string b = 6162 ... 6n be defined 
as 

<Z/b(p,g,r) (g) 

where is the number of occurrences of symbol i in b. If A is a set of strings, 

we define density polynomial associated with A as 

'l'A{p,q,r) = '^'F^{p,q,r). (7) 

aeA 


One can easily show (in a manner similar as done in 0 ) that if one starts 
with a bi-infinite string of symbols drawn from Bernoulli distribution where 
probabilities of 0,1 and 2 are, respectively, p, q and r, then the proportion of sites 
in state k after n iterations of rule / is given by (p, r). This quantity 

will be called density of symbols k after n iterations of /. 

For 3-state rule 140 defined by eqs. (|^ and Q, we generated sets of n- 
step preimages of 0, 1 and 2 for n varying from 1 to 7. Using AT&T FSM 
Library [1] , we constructed minimal finite state machines (FSM) generating these 
sets, and we found that these FSM exhibit regularities which can be exploited 
to produce general expressions for •Ff-,x(j,)(p, g, r). Results are described below. 
Formal proofs are omitted for lack of space, but they are available upon request 
and will be published elsewhere. 


2 Structure of preimages of 1 

The set of n-step perimages of 1 can be described by a finite state machine 
(FSM) schematically shown in Figure The FSM has four parts, denoted by 
P,Q,R and S. Parts Q and S are always the same, while parts P and R consists 
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Fig. 2. Finite state machine for for 6-step (a) and (b) n-step preimages of 1 under the 
rule 140. 


of repeated graph fragments, where the number of repetitions is, respectively, 
n — 2 and n — A. 

Let us consider part Q first. If we start from the leftmost node of P and 
wish to end at the node labeled by a, the only path with such property is 222, 
corresponding to density polynomial . Similarly, if we want to end at &, the 
possible paths are *11 and 221, yielding density polynomial {p-\- q + r)q^ + r'^q. 
For the node c, the possible paths are 021,121, 001,101 and 201, so the density 
polynomial is 2pqr A- q^r + p^q + q^p. Finally, for node d, the paths are 022 and 
122, so that the polynomial is pr'^ + qr^. Let us, therefore, define a vector with 
entries corresponding to density polynomials of paths ending at a, 6, c, and d. 


Q 


{p + q + r)q‘^ + r'^q 
2 pqr + q^r + p'^q + pq^ 

p^2 _|_ ,..2g 


( 8 ) 
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In a very similar fashion, we can construct a vector holding density polynomials 
for paths starting at nodes of segment S labeled, respectively, a, 6, c, and d, and 
ending at the rightmost node. 


5 = 


+ q'^r + 2 {p + q + r) qr + r'^q 
q^ + q^r + [p + qr) qr + r {p + qr)^ 
{p + q + rf 

r'^q -\- {p + q + r)qr -\- {p + q + r)^ q 


(9) 


It is easy to verify that the above components of S correspond to paths starting 
from a (III, 112,12*, 21*, 221), from b (111, 112,12*, 2 * *), from c (***) and 
from d (221, 21*, 1 * *). 

Let us now analyze segment R of the FSM shown in Figure First let us 
suppose that there is no repeated part in segment R, as it would be for the case 
of n = 4, when the number of repetitions is n — 4 = 0. Let us construct a matrix 
R such that Rij represents the density polynomial of all paths starting from 
node j of segment Q and ending in node i of segment S, where i,j G {a, &, c, d}. 
This matrix has the form 


R = 


r 0 0 0 

qq 0 0 

Orp+q+rq 


0 0 0 r 


( 10 ) 


As we can see, the only non-zero entries are diagonal ones and Rh^a = q, Rc,b = r, 
Rc,d = d- This is because there are only three ways to finish at a different node 
that we started, namely if we start from a and finish at b (generating symbol 1 
along the way), if we we start from b and finish at c (generating 2), or if we start 
from d and finish at c (generating 1). 

Suppose now that the repeated fragment in segment R is repeated n — 4 
times. It is not hard to see that density polynomials for all paths from node j 
of segment Q to node i of segment S will be represented by entries of matrix 
jin-3 Furthermore, all paths from the beginning of segment Q to the end of 
segment S will be represented by the density polynomial given by S^R'^~^Q, 
where T denotes transposition (row vector). Since the segment P is represented 
by (p + <7 + the final expression for the density polynomial of n-step 

preimages of 1 is 


S'f-d) (P, g, r) = (p + q + r)"-25^R"-3Q. (11) 
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In order to obtain more explicit expression for q, r), we will need to compute 
jin-3 r ^ q, R is invertible, and one can diagonalize it, 


where 


R = L 
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This yields, after simplification, 


( 12 ) 


(13) 


pq^ {-pr + pq ^ q^) {q\Y 

= —A2(p+ ,.)(,-,■)- 

qr i^—p^r + p^g + pq^ — 2 pqr + — g^r) (rA)" 

A2 (p + g) (g - r) 

g (p^ + p^g + 2p^r + pr^ + 3pgr + r^ + r^g + g^r) A^ 

A (p + r) (p + g) 

where we used A = p + g + r. 

When r = g, matrix R becomes singular. In can be written in Jordan form 
as 

R = L{M + N)L-\ 

where 
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(16) 

Matrices M and N commute, and matrix N is nilpotent, iV^ = 0. Because of 
this, for any integer k, 
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and finally 


'Pf^Hi){p,q,q) = {p + q + qT-^S^ = (^ + 29)-^^^ L(M + 7V)"-"L-^Q 

■(p + 2g)””^ 0 0 o' 

0 q”-2 g"-3 (n - 2) 0 

0 0 g”-2 0 

0 0 -”-2 


= {p + 2qr-^Sn 


0 


q 


L-^Q (18) 


After simplification this yields 


pq^ {n + 1) (qX)'^ q'^ {2p^ + Ap'^q + pq'^ - 2q^) {qXY 


+ 


(q + p) A2 

(p3 + + 4p<7^ + 3 q^) qX^ 


X{q+py 


(19) 


where, as before, X = p + q + r = p A- 2q. 


We shall add here that even though eq. (141 and (191 were derived assuming 


n > 4, they happen to be correct for n = 1, 2 and 3 as well. Let us also remark 
that when one substitutes p = 1, q = 1 and r = 1, then •Z'f-,.(i)(l, 1,1) counts the 
number of preimages of 1. This yields a sequence exhibiting linear-exponential 
growth, 

^ ^ 7 \ 11 

( 20 ) 


WM,l)=()( + L^3” + li9” 


3 Structure of preimages of 0 and 2 


For preimages of 0, FSM generating preimage sets are quite similar as for preim¬ 
ages of 1, thus we will omit details. Similar analysis as in the previous section 
yields, for r g, 

(—pr+pq + q^)pq‘^{qXY pr (—r'^p + q^p + q^ — — qr'^) (rXY 

= —x“-“TP-',)—+-- 

{^p^ + 2p'^q + 2p'^r + 2 r^p + 3 qpr + 2 q^p -f -|- 2 qr^ -|- + q'^r'^ 


and for r = q, 


^f-^(o){P,q,q) = 


(p -f g) (p -f r) A 


(p3 -I- 4p^g -I- 7 q^p + 5 g3) pq^ (n + 1) (gA)*^ 


( 21 ) 


X{p + qf A2 (p -F g) 

q^p (3p2+8pg + 6g2) (gA)' 


(p-fg) A2 


. ( 22 ) 
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Fig. 3. Finite state machines for 4-step (a) and n-step (b) preimages of 2 under the 
rule 140. 


Again, when p = q = r = 1, the density polynomial counts preimages of 0, and 
we obtain 

+ (23) 

This sequence, similarly as the number of preimages of 1, exhibits linear-exponential 
growth. 

For preimages of 2, preimage sets have much simpler structure, shown in 
Figure Density polynomials for them are given by 

= (g+p)rA2”"^ + y (t-A)”. (24) 

The number of preimages, obtained by taking p = g = r = 1, is in this case 

= + (25) 

thus no linearity is present. 

4 Conclusions 

As mentioned in the introduction, density polynomials (p, g, r) represent 

probability of occurence of k after n iterations starting from a Bernoulli distri¬ 
bution with probabilities of 0, 1 and 2 equal to, respectively, p, g, and r, where 
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p + (7 + r = l. Ifwe start with a symmetric Bernoulli distribution where r = q, 
the probability of occurence of 1 after n steps, to be denoted by P„(l), will be 
given by eq. (19) in which we substitute r = q and g = (1 — p)/2. This yields, 
after simplification, 


P„(l) = Poo(l) - ^ {p^n-p^ -pn- 5p- 3/ + l) , 

4(l+p) V ^ / 


where 


Poo(l) = 


(1 - p) (p^ + 5p^ - p + 3) 
4(1 +p)^ 


(26) 


(27) 


It is clear that for 0 < p < 1, Pn(l) tends to Poo(l) as n —>■ oo, and that 
the convergence is linear-exponential in n. Such “degenerate” convergence takes 


place for probability of occurence of 0 as well, as seen in eq. (22). When r ^ q 


in the initial Bernoulli distribution, the convergence is purely exponential, as in 
eq. (fldl) and (pTl). Probability of occurence of 2 is also always exponential, and 


degeneracy is not possible in this case. 

The example of 3-state rule presented here is an interesting instance of a 
phenomenon similar to degenerate hyperbolicity in finite-dimensional dynamical 
systems. It is hoped that it stimulates further research on hyperbolicity in CA. 
The method of constructing density polynomials by using finite state machines 
appears to be quite fruitful, and it should be applicable to many other cellular 
automata rules. 
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